Time-of-flight patterns of ultra-cold bosons in optical lattices 
in various Abelian artificial magnetic field gauges 
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We calculate the time-of-flight patterns of strongly interacting bosons confined in two-dimensional 
square lattice in the presence of an artificial magnetic field using quantum rotor model that is inher- 
ently combined with the Bogolyubov approach. We consider various geometries of the magnetic flux, 
which are expected to be realizable, or have already been implemented in experimental settings. The 
flexibility of the method let us to study cases of the artificial magnetic field being uniform, staggered 
or forming a checkerboard configuration. Effects of additional temporal modulation of the optical 
potential that results from application of Raman lasers driving particle transitions between lattice 
sites are also included. The presented time-of-flight patterns may serve as a verification of chosen 
gauge in experiments, but also provide important hints on unconventional, non-zero momentum con- 
densates, or possibility of observing graphene-like physics resulting from occurrence of Dirac cones 
in artificial magnetic fields in systems of ultra-cold bosons in optical lattices. Also, we elucidate 
on differences between effects of magnetic field in solids and the artificial magnetic field in optical 
lattices, which can be controlled on much higher level leading to effects not possible in condensed 
matter physics. 

PACS numbers: 67.85. Hj, 03.75.K.k, 05.30.Jp 



I. INTRODUCTION 

The experimental observation of Bose Einstein conden- 
sation [T] of trapped atomic gases catalyzed a very large 
further research activity in studies of behavior of atoms 
obeying Bose-Einstein statistics. This led to loading of 
ultra-cold bosonic atoms into optical lattices, which offer 
a clean setting for quantitative and highly precise inves- 
tigations of quantum phase transitions in the strongly in- 
teracting atomic systems [2] . The behavior of the atoms 
bears resemblance to the physics of strongly-correlated 
electronic systems like high-T^ superconductors, which 
are described by similar microscopic Hamiltonians (Hub- 
bard model) [3]. Although the particles that are loaded 
to optical lattices are electrically neutral, it is possible to 
impose additional external potential, which forces them 
to behave exactly like charged particles interacting with 
an external magnetic field. In the simplest case the po- 
tential can result from rotation, following from formal 
equivalence between the Lorentz force and the Coriolis 
force [HIS]. However, more control is obtained using ad- 
ditional photon-assisted tunneling to coherently transfer 
atoms from one internal state to another. This induces 
a non- vanishing phase of particles moving along a closed 
path, which simulates magnetic flux through the lattice 
[SHS]. Such techniques are under very active investiga- 
tion. Since, quantum optics technology provides unprece- 
dented degree of manipulation of structure of such im- 
posed magnetic flux, it allows for obtaining very strong 
magnetic flelds for neutral atoms: both Abelian [9, ID] 
and non- Abelian |1H [12]. A long-term goal of these 
experiments is to achieve the quantum Hall regime, in 



which very high value of effective flux opens an avenue 
to study the effects that are not achievable in the con- 
ventional solid state physics [151 [T^ . 

The dynamic properties of the trapped ultra-cold 
atoms can be investigated in time-of-flight (TOF) ex- 
periments, in which the trapping potential is suddenly 
switched off. As the atoms are no longer being localized 
spatially, their scattering ratio decreases and the momen- 
tum distribution becomes temporarily frozen, with grav- 
ity being the only force acting on the atomic cloud. The 
infra-red absorption images taken after arbitrary expan- 
sion time show locations of the atoms, which are directly 
related to the the distribution of the momenta n (k) in 
the system before the potential was switched off. In the 
superfluid state (SF) the TOF images exhibit character- 
istic sharp maxima related to long-range phase coherence 
of the condensate. While the phase fluctuations are be- 
ing increased by stronger interactions between atoms, the 
sharp features disappear and n(k) becomes a wide maxi- 
mum indicating the presence of the Mott insulating (MI) 
state [2,. 

Theoretical challenge in describing the TOF patterns 
results from the dynamic nature of the problem: deter- 
mination of the momentum distribution n (k) requires 
precise knowledge of spatial correlations between atoms, 
namely atom-atom correlation function. This precludes 
use of methods based on mean-field approximation. In 
the present paper we apply a recently proposed combi- 
nation of quantum rotor (QR) approach and Bogolyubov 
method, which has been successfully applied to investi- 
gate correlations in systems of cold atoms in optical lat- 
tices (e.g. time-of-flight patterns |T5], spectral functions 



|16j). It is also a natural extension of the QR model that 
was used to describe the phase diagram, also in the pres- 
ence of the artificial magnetic field llTffTO] . The QR ap- 
proach has been verified pi^.f^ri using other methods, like 
Monte Carlo numerical calculations [22 and diagram- 
matic perturbation theory |23| . The QR phase diagrams 
were also analyzed in the context of an analytical works: 
mean-field theory [24 and Fade analysis [25,. 

It is our goal to calculate the TOF patterns for vari- 
ous gauges of the magnetic field that are expected to be 
realizable or have already been implemented in the exper- 
imental settings to investigate the dynamics of the con- 
densate and the phase transition to the localized state. 
The remainder of the paper is as follows: in. Sec II 
we introduce the model Hamiltonian relevant to strongly 
correlated bosons confined in two-dimensional square lat- 
tice. Then we apply the synthetic magnetic field, which 
modifies the hopping term in the Hamiltonian and add 
additional temporal modulation of the optical potential. 
The main points of our approach that lead to calculation 
of the atom-atom correlation function are summarized 
in Sec. III. Furthermore, we present the time-of-flight 
patterns in the following Section. Our results are sum- 
marized in Sec. V, while the dispersion relations and the 
resulting lattice densities of states used in calculations 
are presented in Appendix. 



II. MODEL HAMILTONIAN 

The essential physics of bosons in optical lattice can be 
captured using the single-band Bose-Hubbard model. In 
this description, the particles move within a tight-binding 
scheme and interact only through on-site repulsion re- 
sulting from interatomic collisions (since the atoms are 
neutral). The Hamiltonian is given by. 
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where bi {b\\ is the boson destruction (creation) opera- 
tor at a site «, ni = b\bi is the density operator, C/ > 
is the on-site repulsion and /i is the chemical potential, 
which controls the number of bosons. Here, (j, j) denotes 
summation over the nearest-neighbor sites. Finally, tij is 
the hopping matrix element, which is non-zero only for 
the nearest neighbors and equal to t. For any given lat- 
tice geometry and depth, both t and U can be calculated 
directly by finding the respective Wannier function basis 
[26] . Introduction of the synthetic magnetic field B (po- 
tential that acts on neutral particles in the same fashion 
as the magnetic field acts on charges, e.g. rotation of the 
system, laser stirring, selective driving of hopping with 
Raman lasers) leads to introduction of the Peierls phase 
factor: 



exp 



27rj 

$0 



A-dl 



(2) 



which is a consequence of the gauge invariance of the 
Schrodinger equation, where B = V x A (r) and <i>o = 
h/e is the flux quantum, with A (r) being the vector 
potential. This leads to modification of the hopping term: 
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that as a result, instead of being a real value, becomes a 
complex number: 
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Furthermore, particles hopping along closed loops of the 
lattice cell (area a) gain an additional phase (/) = 27r/ im- 
posed by the external uniform synthetic magnetic field 
potential, where / = aBe/2TTh. We also permit for 
additional spatial modification of the on-site potential 
( — 1)' A, which allows us to describe the effect of tem- 
poral modulation by the photon-assisted tunneling that 
are used to drive the phase change in some experiments. 
The range of actual gauges (shapes of the vector poten- 
tial A) applied to the system, which can be realized ex- 
perimentally, is very wide (see. Sec. IV). The change 
of the hopping parameter in Eq. Q also modifies the 
band structure, which becomes very complex. Compli- 
cated multi-band dispersion relations provide difficulties 
in calculating analytical formulas in the uniform case for 
the lattice density of states (DOS) limiting availability 
to a few selected values of / |17| . 



III. CORRELATION FUNCTIONS 

In optical lattices, the phase transition between su- 
perffuid and Mott insulator states occurs in the regime 
of intermediate to strong interactions {U ^ t). As a 
result, a theory that goes beyond standard Bogoliubov 
approximation is required. To this end, we calculate the 
one-particle correlation function that is necessary to pre- 
dict the time-of-flight patterns using the quantum rotor 
approach (see, Ref. [19]) combined with the Bogolyubov 
method that has been recently proposed and succesfuUy 
applied to systems of bosons in optical lattices [15] . This 
scenario provides a picture of quasiparticles and energy 
excitations in the strong interaction limit, where the tran- 
sition between the superfluid and the Mott state is driven 
by phase fluctuations. The approach is based on sep- 
aration of the problem into the amplitude of the Bose 
field and the fiuctuating phase that was absent in the 
original Bogoliubov problem. As a results, one arrives 
at a formalism, where the one-particle correlation func- 
tions are treated self-consistently and permit us to inves- 
tigate a whole range of phenomena described by the Bose- 
Hubbard Hamiltonian. Furthermore, the phase fiuctua- 
tions are described within the quantum spherical model 



|27| . which goes beyond mean- field approximation includ- 
ing both quantum and spatial correlations. Although, 
the approach easily allows for non-zero temperatures, in 
the following we restrict ourselves to the description of 
the ground state of the system (T = 0) . As the details of 
calculations have been extensively presented in Ref. |15j . 
we only summarize the main steps of the approach here. 
We start by introducing the functional integral represen- 
tation of the model in Eq. (fTl) in terms of the the complex 
fields Oi (r), which leads to the partition function: 
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with the action S given by 
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where /3 = l/k^T and T being temperature. Next, 
we perform the local gauge transformation to the new 
bosonic variables 



a* (r) = h (r) exp [iip^ (r)] 



(7) 



It allows to extract phase variable ipi (r), which order- 
ing naturally describes the superfluid - Mott insulator 
transition, and the amplitude bi (r) that is related to the 
superfluid density. As a result, the the partition function 
becomes: 
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with the action S \b, b,(p\ =5' [a, a] . The statistical sum 
in Eq. dsl can be integrated over the phase or amplitude 
variables leading to phase-only of amplitude-only actions: 
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The main point of the approach is the calculation of the 
action S^p [(p] in Eq. dol, which describes the phase-only 
model with amplitudes integrated out. It is subsequently 
mapped onto the quantum spherical model, which can be 
solved analytically. 

As a result of the variable transformation in Eq. 
(l7|, the superfluid order parameter, which non- vanishing 
value signals a macroscopic quantum phase coherence 
(identifled as the superfluid state), factorizes: 
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This reflects the fact that all atoms in the condensate 
form a coherent matter wave having the same phase. The 



averages in Eq. (11) are deflned as: 
/ , _J[Dx]...e- 



IP 



x\ e 



-5,N 



(12) 



for X ~ a, b, (f and the respective actions: S [a, a], S [b, b\ 
or S^ [ip]. Furthermore, we parametrize the boson flelds 



b^ (r) = &o + 6. {t) , 
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where bg = \/No is the Bose condensate macroscopic oc- 
cupation and &^ (r) is the amplitude fluctuation around 
the mean value bo. As a result, the superfluid order pa- 
rameter becomes: 



*B = boniQ, 
where tuq is phase order parameter: 

mo = (exp [iipi (r)]) 
The atom-atom correlation function 



(14) 
(15) 



C.,{T) = {a,{T)a,{T))^ (16) 

also factorizes due to the variable transformation in Eq. 
(I?! becoming: 

C^J (r) = (b'b,)^ (exp [ip, (r) - p, (r)])^ (17) 

with the averages that can be calculated analytically for 
any lattice, for which the dispersion relation ik (Fourier 
transform of the hopping tij) is known [15J. The momen- 
tum distribution of the atoms in optical lattice is then a 
Fourier transform of the correlation function: 
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This leads to the density of particles in the time-of-flight 
experiments |15l [5S] : 



n(r) = 



ht, 



T4^ k = 



ht„ 



k = 



htp 



(19) 

where \W (k)| is the envelope of the Fourier transform of 
the Wannier function for the chosen optical lattice and 
te is the expansion time. It should be pointed out that 
the envelope |Ty(k)| can in principle depend not only 
on the optical lattice potential, but also on the presence 
of the artiflcial magnetic fleld [29_. However, since its 
calculation goes beyond the scope of the present work, 
we use the standard form: 
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where Wt = M^/mwo with wq being the size of the on- 
site Wannier function [T5] . This choice can be justified by 
comparing the resulting TOF patterns with experimental 
ones [lOj and observing the conformity of the particle 
density decays as a function of r in both cases (see. Sec. 
iVDi. 



a} uniform magnetic field b) staggered potential 



c) checkerboard staggered field d) uniaxially staggered field 
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Figure 1: (Color online) Geometry of the artificial magnetic field in various gauges resulting in different flux configuration: 
a) uniform, b) uniform with additional staggered potential, c) staggered flux with checkerboard arrangement, d) uniaxially 
staggered flux. Value of the flux per plaquette being the phase acquired by a particle traveling around an elementary cell is (j). 
Additional staggering on-site potential is denoted by ±A. 



IV. EXPERIMENTAL GAUGES OF THE 
SYNTHETIC MAGNETIC FIELD 

The monicntum distribution is an important observ- 
able since it allows to identify whether the atoms in the 
optical lattice are in the superfluid or Mott insulating 
state. Although, the existence of sharp peaks in the time- 
of-flight images that used to be unequivocally associated 
with the emergence of the superfluidity is not believed 
to be sufficient criterion [23 , recent analysis have shown 
that the momentum distribution can be used to make 
pretty accurate estimations about location of the critical 
regime [50| . 

In the following, we calculate the time-of-flight pat- 
terns resulting from Eq. (19 I for various gauges of the 
artificial magnetic field. We start from determining the 



dispersion relation t\^ for the chosen gauge, which allows 
us to obtain the atom-atom correlation function in Eq. 
( 17 1 by using the procedure described in details in Ref. 

We present all the results along similar scheme: first 
we plot tk for the chosen gauge and then the calculated 
time-of-fiight patterns in the superfluid state, near the 
SF-MI phase transition and in the Mott insulator. We 
measure the interaction strength between atoms using 
the experimental quantities Vq and Er instead of t and 
U , where Vg is the optical potential depth and E]:j - the 
recoil energy. The relation of Vo/En to t/U is presented 
in Ref. [T5] . One should also note that the minimum 
of kinetic energy in the Hamiltonian in Eq. (ITl corre- 
sponds to the maximum of it due to the minus sign in 
the hopping term. 



A. Uniform magnetic field 



We start with the uniform artificial magnetic field, 
which acts on atoms in optical lattices in identical way 
as a homogeneous magnetic fields acts on electrons in 
solids (see. Fig. flk). Every elementary cell of the lat- 
tice is pierced by a fraction / of the elementary fiux, 
which leads to additional phase (j) acquired by a parti- 
cle moving around the cell equal to cf) = 27r/. Such a 
configuration of the artificial magnetic field can be real- 
ized using various gauges, e.g. Landau A = i3(0,j/,0), 
or symmetric A = ^ (— y,a;,0). This results in increase 
of the elementary cell, since translational symmetry is 



locally broken for non-integer values of /. If / is a ra- 
tional, being equal to / = p/q, the cell enlarges q-fold, 
while the Brillouin zones shrinks by the same factor (lead- 
ing to so-called magnetic Brillouin zone). The quasipar- 
ticle spectrum has a complicated multi-band structure 
known as the Hofstadter butterfiy [14J (the denomina- 
tor q determines the number of sub-bands) and can be 
generated using Harper's equation fT3. Although, the 
general solution is unknown, for special values of / equal 
to / = 1/2, 1/3, 1/4, 1/6, 1/8, 3/8 both dispersion rela- 
tion and lattice density of states have been analytically 
calculated jH]. We see, that the denominator of the ex- 
pression describing the magnetic field / = p/q determines 
the number of bands. 
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Figure 2: (Color online) Dispersion relation and time-of-flight 
patterns for atoms in optical lattice without (a) and under 
uniform artificial magnetic field: b) / = f/2, c) / = f/3, 
d) / = 1/4. The TOF images depict superfluid state (the 
first column), system close to the phase transition (the second 
column) and the Mott insulator (the last column). 



The dispersion relations for / = 0, 1/2, 1/3 and 1/4 
have been presented in Fig. |2] along with the resulting 
TOF patterns. In the superfluid state, the spectrum has 
sharp coherence peaks, which slowly fade away, when the 
system is driven towards the phase transition. Finally, in 
the Mott insulator, the TOF patterns become a wide, 
feature-less maximum with no signatures of the phase 
coherence. The most notable influence of the synthetic 
magnetic field is a change in periodicity leading to shrink- 
ing of the Brillouin zone and accompanying dense packing 
of the coherence peaks in the SF state. As a result, small 
fluxes could be difficult to detect. In the typical TOF ex- 
periment the 2'hk of the Brillouin zone corresponds to the 
50 pixels on the charge-coupled-device (CCD) camera. 
The smallest detectable separation between two momen- 
tum peaks is around 10 pixels which matches the 1/5 of 
the Brillouin zone, therefore the fluxes below (j) < 2tt/5 
(/ = 1/5) will be hardly recognizable from the experi- 
mental data. Moreover, the key role is the proper prepa- 
ration of the ground state, in which the coherence over 
large area of the real space is obtained to avoid further 
peaks broadening [31J. In the theoretical calculations the 
resolution of the TOF diagrams can be in principle arbi- 
trarily high and two peaks will become indistinguishable 
when the distance between them becomes of the order of 
the full width at half maximum of the n (fc) peaks, which 
in this case is about 3.5% of the flrst Brillouin zone width. 



1. f = l/2(<l, = n) 

Since, the single-particle spectrum is symmetric 
around / = 1/2: tk(/) = ^k(l ^ /)■ the strongest possi- 
ble uniform artiflcial magnetic fleld that can be achieved 
is / = 1/2, which results in flux (j) — t: per plaque- 
tte. In this case ik has two sub-bands, which meet at 
tko — forming Dirac cones (the spectrum is linear near 
li-D and rotationally symmetric: t^ ^ |k — k/jj). This 
leads to potential possibility of observing graphene-like 
physics in optical lattices. The time-of-flight patterns 
for / = 1/2 are presented in Fig. l3] In the case when 
bosons are free to occupy any of the sub-bands (see. Fig. 
pa) , the TOF images show coherence peaks in the super- 
fluid state at the momenta, for which the kinetic energy 
assumes minimal values (k = {nn, mn} , where n, m are 
integers, which are basically k = point repeated by 
periodicity of the reciprocal lattice) . The k/j points cor- 
responding to intersections of the bands do not show in 
the TOF patterns. However, if all particles were occu- 
pying only the upper band, the kinetic energy minima 
would appear at k^i, which has been presented in Fig. 
|3]3. The resulting superfluid state exhibits non-zero mo- 
mentum, with k = component totally removed. On 
the other hand, in case of populating the lower band 
only (see. Fig. l3t), the resulting picture does not differ 
much from the scenario, when occupation of both bands 
is allowed (in Fig. [3k), however some slight differences 
are noticeable. It results from the fact that although we 
are investigating the ground state of the system, not all 
bosons occupy the lowest energy state. Since the parti- 
cles are interacting, only a fraction of them contributes to 
the condensate (thus occupies the lowest energy state), 
while the rest can be driven to higher energy states by 
the quantum ffuctuations, which are present even in zero 
temperature. This also leads to a conclusion that in sys- 
tems of bosons the properties of the superffuid state are 
determined by the points in the fc-space around the min- 
ima of the kinetic energy. However, in the vicinity of 
the SF-MI phase transition the sharp maxima are gone, 
and the TOF images depict momentum distribution of 
the incoherent particles. Surprisingly, they contain weak 
maxima around k = points regardless of location of the 
superffuid peaks. 

In order to observe effects resulting from existence of 
the Dirac cones in the excitation spectra, population of 
the respective bands has to be engineered. It can be ex- 
perimentally realized using projection of the condensates 
onto a desired Bloch state [ST: the system of bosons, 
which initially is in the superffuid phase, is released from 
the trap and expands freely for a short period of time. 
Then, a moving optical lattice is introduced, which is cre- 
ated by laser beams with additional acousto-optic modu- 
lators that allow for shifting positions of the lattice min- 
ima. As a result, depending on the modulation, the BEG 
can be loaded to a lattice state with an arbitrary and 
well-deffned quasi-momentum. Using this approach, it 
was possible to access different energy bands of the *^Rb 
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Figure 3: (Color online) Dispersion relation and time-of-flight 
patterns for atoms in optical lattice under uniform artificial 
magnetic field / = f/2: a) both sub-bands filled, b) higher 
sub-band filled, c) lower sub-band filled. The TOP images 
depict superfluid state (the first column, Vo/Er = 5), system 
close to the phase transition (the second column, Vo/ Er = 
7.5) and the Mott insulator (the last column, Vo/Er — f5). 



atoms allowing a high precision studies of the lensing ef- 
fect on a Bose-Einstein condensate |33| . 

It should be also noted that the appearance of the 
Dirac intersections in the Hofstadter spectra occurs for 
magnetic fields / — p/q, for which q is even and is never 
observed for odd values Fig. [2] 



B. Uniaxial staggered potential 
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Figure 4: (Color online) Dispersion relation and time-of-flight 
patterns for atoms in optical lattice under uniform artificial 
magnetic field / = 1/2 and additional staggering potential A 
driving particle hopping for A equal to: a) 0, b) 0.1, c) 0.6, d) 
1.2, e) 2. The TOF images depict superfluid state (the first 
column, Vq/Er = 5), system close to the phase transition 
(the second column, Vq/Er = 7.5) and the Mott insulator 
(the last column, Vq/Er = 10). 



Application of the additional staggering potential 
A/2t = A that drives hopping of atoms between cho- 
sen lattice sites [see Sec. II] is a natural extension of the 
system in the uniform magnetic fields (see, Fig. ma). It 
allows to manipulate the Dirac cones: change the dis- 
tance between them in the fc-space and annihilate them 
when two of them merge [33]. For / = 1/2 (flux through 
the elementary cell (p = ir) and A = 0, the system is iden- 
tical to described in Sec. |IV A lj _and the resulting TOF 
patterns are presented in Fig. |4^. While the A is be- 
ing increased, the Dirac points move closer to each other 
(see, Fig. kp-c) and for A = 1 - merge. For A > 0, the 
Dirac points annihilate and the single-particle spectrum 
becomes gapped (Fig. pH-e). However, as in the previous 
case of the uniform field, if both sub-bands are populated, 
the condensation of bosons occurs around the bottom of 
the lower band, thus Dirac cones have no effect on the 
TOF patterns. On the other hand, increasing value of 
A strongly enhances hopping along one direction leading 
to slow decline in weight of k = maximum, enlarging 
(uTT, ±7r) components (with n ^ 0). This effect can be re- 
versed by enlarging interatomic interactions (U/t): near 



the superfiuid - Mott insulator phase transition the mo- 
bility of atoms is naturally decreased. Once again, this 
allows to observe the weak maxima in momentum distri- 
bution of incoherent particles which are located around 
k = regardless of the position of the superfiuid phase 
coherence peaks. 

A slightly different behavior can be observed for a sys- 
tem without artificial magnetic field, but with the stag- 
gering potential A (see. Fig. p|. Since the fiux is missing, 
the maxima occur at k = 0, however their density in k^ 
direction is doubled due to increased size of an elemen- 
tary cell (with the width of the Brillouin zone halved). 
Increase of A leads to smearing of the peaks in k^ di- 
rection. The dispersion relations are quite different: the 
Dirac cones are not present and the spectrum is gapped 
for every value of A. Also, a noticeable difference occurs 
for A « 0: in the system with / = the TOF images 
change discontinuously while going from non-zero to zero 
value of A. This results from the fact that presence of 
the staggering potential breaks the translational symme- 
try doubling the size of the elementary cell for every value 
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Figure 5: (Color online) Dispersion relation and time-of-flight 
patterns for atoms in optical lattice under additional stagger- 
ing potential A driving particle hopping for A equal to: a) 0, 
b) 0.6, c) 1.2, d) 2. The TOF images depict superfluid state 
(the first column), system close to the phase transition (the 
second column) and the Mott insulator (the last column). 



of A, but not for A = (see, differences between Figs. |5| 
and b-e). This is in contrast to the / — 1/2 case, where 
the enlargement of the elementary cell resulting from the 
presence of the tt flux per plaquette and the staggering 
potential are the same. 
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Figure 6: (Color online) Dispersion relation and time-of-flight 
patterns for atoms in optical lattice in checkerboard staggered 
artificial magnetic field with fiux (f> per plaquette: a) 0, b) 
TT, c) 27r. The TOF images depict superfluid state (the flrst 
column), system close to the phase transition (the second col- 
umn) and the Mott insulator (the last column). 



Wp) ) the energy minima of the two conical points be- 
come equal, leading to peaks in momentum distribution 
located in k = and k = (±7r, ±7r) points. Finally Fig. 
l6l;) , for larger fluxes the non-zero momentum state takes 
over reaching maximum intensity for (p = 2-k. It should 
be stressed that in a naive view the impact of the flux 7i27r 
(n being integer) should be negligible. However, here the 
(j) = ±27r flux leads to non-trivial superfluid phase with 
non-zero momentum. This issue will be discussed in the 
following subsections. 



C. Checkerboard staggered flux 



D. Uniaxially staggered flux 



A time-independent lattice model with an artiflcial 
staggered magnetic field that is used in the present 
work can effectively describe time-dependent optical lat- 
tice with staggered particle current in the tight-binding 
regime [36]. As a result, it is possible to describe group 
of experiments that use temporal modification of the 
optical potential. They allow reaching regimes, where 
anisotropic Dirac cones emerge in the single-particle spec- 
trum leading to two inequivalent conical points in the en- 
ergy band, which results in two distinct energy minima 
that depend on the magnitude of the staggered magnetic 
fiux (j). Consequently, it is possible to realize the artifi- 
cial magnetic field as presented in Fig. [l]:, where the fiux 
is staggered and arranged in checkerboard configuration. 
This method allows to reach values of fiux per plaque- 
tte ranging from — 27r to 27r. The resulting time-of-fiight 
images are presented in Fig. |6] For small fiuxes </> the 
effect of the magnetic field is hardly noticeable, as the 
strong maximum at k = is visible. For (f) = tt (see Fig. 



Photon-assisted tunneling in an optical superlattice 
generating large tunable effective magnetic fields for 
ultra-cold atoms demonstrated possibility of realization 
of the large tunable uniaxially staggered field (where the 
spatial average of the fiux is zero) ^lOj . It was shown that 
the atomic sample relaxes to the minima of the magnetic 
band structure, realizing an analogue of a frustrated clas- 
sical spin system. The obtained time-of-fiight patterns 
for various system hopping anisotropics |10[ [35] agree 
well with the ones calculated with the method presented 
in the current work (see. Fig. ul. Positions of the max- 
ima of the momentum distribution are correctly recreated 
as well as the decay of the envelope of the TOF image 
substantiating choice of the module of the Wannier func- 
tion in Eq. (20). For isotropic system, the time-of-fiight 



patterns exhibit two minima located around the zero mo- 
mentum at k = ± (7r/4, 7r/4). While the anisotropy be- 
tween raw hopping in the x (i^+ij) and y (tij^i) di- 



rection is introduced, for i. 



/U 



?+lj / ^i.j-\-l 



< \/2, the peaks 
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Figure 7: (Color online) Dispersion relation and time-of-flight 
patterns for atoms in optical lattice in uniaxially staggered 
artificial magnetic field with flux per plaquette (j> = n/2 for 
hopping anisotropy leading ratio of fi+ij/ii,j+i equal to : a) 
1, b) 2, c) 3. The TOF images depict superfluid state (the 
first column), system close to the phase transition (the second 
column) and the Mott insulator (the last column). 



split into pairs of peaks in agreement with changes of the 
magnetic band structure. 



E. Arbitrary gauge geometry 

Although, the flux configuration in the case of the uni- 

and Fig. IT]) is 
7r/2 and 



IV D 



axially staggered flux (see. Sec 
pretty regular: uniaxially alternating values of 
— 7r/2 every second plaquette, the structure of the TOF 
images is complicated and strongly dependent on the lat- 
tice parameters like hopping anisotropy. It results from 
complex gauge that was used in the experiment in Ref. 
[lOj. However, the same configuration of the flux per 
plaquette can be obtained for much simpler gauge, as 
presented in Fig. ^p (the dispersion relation obtained in 
the same manner as in Ref. [lOj, however without phase 
change along k^ direction hopping). The comparison of 
the resulting time-of-flight patterns is presented in Fig. 
|9] and it is clearly visible that they are not identical. It 
can be seen that in case of atoms moving in the tight- 
binding scheme in the optical lattice the value of the flux 
being assigned to an elementary cell does not determine 
the momentum distribution of the particles. Other than 
that, what is crucial is the change of the quantum phase 
that occurs at every bond that the particle travels along, 
since all the jumps are separated acts rather than a con- 
tinuous move. This can lead to a non-intuitive situation 
of particles exhibiting the influence of the vector poten- 
tial fleld resulting of a non-zero phase change on selected 
bonds, although the total phase change on closed trajec- 
tory around an elementary cell is zero, which also means 
that so is the effective flux per plaquette. Such situation 
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Figure 8: (Color online) Comparison of the gauge used in Ref. 
|10| and a simpler one resulting in the same configuration of 
the fiuxes (denoted by thick, blue arrows located inside the 
unit cell) ±7r/2, while traveling around a plaquette. Thin 
arrows (red) denote direction of hopping along which a par- 
ticle phase changes by n/2 (— 7r/2 in an opposite direction), 
while tubes (green) represent regular hopping with no phase 
change. 



is presented in Fig. 10 the configuration of the gauge 
results in change of phase equal to zero when a parti- 
cle travels around a closed loop. However, the resulting 
time-of-flight patterns are still dependent on the phase 
acquired during a single jump leading to dispersion rela- 
tion: 



t (k) = t [cos {kx 



- cos {ky + 4>)] . 



(21) 



For free condensate (without optical lattice) the super- 
fluid velocity: 



-V0. 



(22) 



The kinetic energy Ek can be expanded around its min- 
imum at k^ = (—0, —0), which leads to Vs ^ dEk/dk — 
Kk^/m. As a result, change of the particles momentum 
resulting from phase acquired on a jump along a single 
bond is simply: 



fc- V(^, 



(23) 



10 It is 



which is consistent with the results in Fig. 
worth to notice that using such simple gauge conflgu- 
ration one obtains a finite momentum superfluid phase 
with the non-zero phase change imposed on the conden- 
sate wave function but with the zero value of the artificial 
magnetic field. 



V. SUMMARY 

In conclusion, we have analyzed the correlations 
between strongly interacting bosons confined in two- 
dimensional square lattice in the presence of an artifi- 
cial magnetic field using quantum rotor model that is in- 
herently combined with the Bogolyubov approach. The 
fiexibility of the method and its sensitivity to the spa- 
tial fiuctuations allows us to consider various geome- 
tries of the magnetic fiux (uniform, checkerboard, uni- 
axially staggered), which are expected to be realizable. 
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Figure 9; (Color online) Comparison dispersion relation and 
time-of-fliglit patterns for atoms in optical lattice in differ- 
ent gauges (see, Fig. [8|, which lead to the same uniaxially 
staggered configuration of the magnetic flux. 
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tides on a single hop is marginal. However, in optical 
lattices, the attainable values of the flux are very high 
and the phase change on single bonds can be individ- 
ually controlled. This leads to the strong dependence 
of the time-of-flight patterns and the atom dynamics on 
the specific gauge configuration rather then the resulting 
magnetic flux. Therefore, in systems of strongly interact- 
ing bosons, only the minima of the kinetic energy (max- 
ima of the dispersion relation in the fc-space) determine 
the superfluid properties of the ultra-cold atoms confined 
in optical lattices. At the same time, the interactions be- 
tween atoms change the phase stiffness and the density 
of the condensate rather than the dynamical properties 
of the coherent particles. 
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Figure 10: (Color online) Gauge resulting in zero flux (red ar- 
rows denote direction of hopping along which a particle phase 
changes by (p) and resulting time-of-flight patterns for atoms 
in optical lattice for Vo/Er — 7.5 and a) (j) — 0, h), (j) = n/3, 
c) <p = 2tt/3 and d) = vr. Dotted boxes mark the first 
Brillouin zone. 



or have already been implemented in experimental set- 
tings. Furthermore, we have calculated the time-of- 
flight patterns, which give information about dynamics 
of condensed atoms surrendered to the artificial magnetic 
field. We validate our approach by successfully recreat- 
ing experimentally observed TOF images and recovering 
the superfluid-Mott insulator phase transition, which is 
driven by the interactions. Also, we show in which con- 
ditions novel superfluid phases with non-zero momentum 
can arise leading to observation of Dirac-like physics in 
optical lattices. Furthermore, we deduce that the crucial 
element for a proper recreation of the time-of flight pat- 
terns is not the flux conflguration but rather the change of 
the quantum phase that occurs at every bond described 
by the dispersion relation and consequently the density 
of states. Thus, the TOF images result directly from 
the chosen and experimentally realized gauge. This is 
in clear contrast with solid state physics, where attain- 
able values of flux per elementary cell are very small (flux 
/ = 1/2 would require flelds of the order of 10^ - lO^T). 
As a result, the spatial change of the vector potential 
(gauge) is very gradual and the phase acquired by par- 



VI. APPENDIX 

Quasiparticle spectrum and the density of states 

In the following, we present dispersion relations and 
the resulting lattice densities of states for various flux 
geometries, which were described in Sec. |IV[ 



1. Staggered potential with uniform flux f=0 for 
two-dimensional square lattice 

The dispersion relation contains two sub-bands P 



to (fc) = 2t (cos ky ± Vcos^fc^ + aA (24) 

and results from the single-particle Schrodinger equation: 



E^„ 



-t^, 



rn,n+l 



tll^v 



ti,, 



rn+l,n 



-iVm-l,« + (-1)" AV'm^n, 



(25) 



which allows to deduce hopping elements tmn {^n, n num- 
ber the lattice sites in x and y directions, E is the energy 
and "0 - the wave function) between neighboring sites. 
The density of states is given by a nonlinear convolution: 

Po {E, A) = ^^-^ I dxpiD {x) 



2ttH 



X PlD 



± Vx^ + A2 



(26) 
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where pijj (x) is a one-dimensional lattice density of 
states: 



PiD {E) 



1 



(27) 



resulting from t{k) = 2tcosfc2: dispersion relation. 



2. Staggered potential with uniform flux f=l/2 for the 
two-dimensional sguare lattice 



The dispersion relation f34J: 



|ii/2 (fc) I = 2iy cos2 k-^ + (cos ky - A) (28) 

for A = is equal to the uniform external magnetic field 
with / = 1/2. It results from the Schrodinger equation: 

- tlprn+l^n - tljjm-l,n + (-1)™ AV'm,„. (29) 

The resulting DOS reads: 



and results from the Schrodinger equation: 



Ei;„ 






-tip, 

tlpm+l,ne 



-le 



^(-1)" 



-«(-iFf_W,.„ ^^P^l-i)"! (32) 



^V^7n—l,n^ 



with p = ni ^ n. The DOS can be written in the form: 



p(i?,0) = 



1^ 

27r2i 



da; 



cos ■ 



-piD [x) 



E^ 



(33) 



X PlD 



where any value of flux (p is allowed. 

4-. Uniaxially staggered flux 



The dispersion relation for gauge configuration used in 
experiments presented in Ref. |10j consists of four sub- 
bands: 



PV.(^.A) -§-Jd. 



X PlD 



PlD (x) 




^{¥y-i--^f 



(x-A)M. (30) 



3. Checkerboard staggered flux for the two-dimensional 
square lattice 

The dispersion relation for the flux (p is given by the 
formula 1361: 



ii,2 (k) = sin fcj; — 77 cos fcy 



± Jt]'^ -2 sin 2kx+r]'^ sin 2ky-^ 2 



t3.4(k) = ^i,2[k-(|,|) 



(34) 



where 77 is a hoping anisotropy ratio, while the dispersion 
used in Fig. [8b, reads: 



ti,2 (k) = cos ky — sin ky — ^2 cos 2kx -I- sin 2ky -f 3 



tzA (k) = ti,2 



k-(|,o 



(35) 



|i0(fc)| = 2< 



2cos(-lcos^^^jcos^^^ 



COS 



2 / '^x ^~ "'y 



COS 



, h — h 
2 / "^a; '''y 



1/2 



In both cases, analytical formulas for the lattice density 
of states cannot be easily obtained. The details of calcu- 
lation of the dispersion in Eq. ( 34 1 can be found in Ref. 



(31) 110|. 
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